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ABSTRACT 

Ixr the study of the propagation of high power microwave pulse, one of the main 
concerns is how to minimize the energy loss of the pulse before reaching the destination. 
In the very high power region, one has to prevent the cutoff reflection caused by the 
excessive ionization in the background air. A frequency auto-conversion process which 
caxx lead to reflectionless propagation of powerful EM pulses in self-generated plasmas is 
studied. The theory shows that under the proper condition the carrier frequency co of the 
pulse will indeed shift upward with the growth of plasma frequency co pe . Thus, the plasma 
during breakdown will always remain transparent to the pulse (i.e., co > co pe )• 

A chamber experiment to demonstrating the frequency auto-conversion during the 
pulse propagation through the self-generated plasma is then conducted in a chamber. The 
detected frequency shift is compared with the theoretical result calculated by using the 
measured electron density distribution along the propagation path of the pulse. Good 
agreement between the theory and the experiment results is obtained. . 

1. INTRODUCTION 

In the applications of high power microwave pulses for radar and directed energy 
system, the effect of air breakdown on the propagation of the pulses through the 
atmosphere has been a central issue of research interest 17 This is because air breakdown 
produces ionization phenomena that can radically modify wave propagation. Moreover, a 
phenomenon called 'tail erosion' plays a primary role in limiting the transmission of 
pulse energy 6 Thus, a considerable effort for such a study has been pursued over the past 
decade. Both fluid and kinetic transport codes have been developed to give a self- 
consistent description of the propagation process. Many laboratory experiments have also 
been conducted to verify the results of computer codes. 

The phenomenon of tail erosion can be described as follows. During the air 
breakdown, a space-time dependent plasma is generated. Considering any point along the 
path of the pulse, the plasma density is building up while the pulse is passing by. Since it 
takes finite time for the plasma to build up, the leading edge of the pulse can only have 
chance to interact with very few plasma electrons and is hardly affected. However, the 
plasma density at that point grows exponentially as the pulse continues passing by. By the 
time as the tail portion of the pulse is about to pass that point, a dense plasma is already 
present and can then attenuate this portion of the pulse significantly. The degree of tail 
erosion usually increases with the intensity of the pulse. By examining the dependence of 
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: -flection signal on the degree of tail erosion in the transmission signal, two mechanism 
responsible for the tail erosion phenomenon have been identified experimentally in our 
, revious work 6 . The first one is the collisional damping by the pulse-generated 
nderdense plasma. It is shown that erosion to the tail of the pulse by this mechanism is 
ot complete and there is no anomalous reflection. The loss of pulse's energy in this case 
;;\ay still be within the tolerable level. The second one is the cutoff reflection caused by the 
; ’ulse-generated overdense plasma screen. In this case the incident pulse has higher 
intensity than that of the previous case and anomalous reflection is detected. The result 
hows that a large portion of the pulse is more or less eroded away during the finite 
propagation period. It is clear that the second mechanism is far more severe, in causing 
energy loss of the pulse than the first one. Once the reflection mechanism is in play, the 
remaining pulse becomes too narrow to deliver enough energy as desired. 

Obviously, one simple way to avoid cutoff reflection is to restrict the power of the 
pulse for a given pulse length to be lower than a critical power level. Unfortunately, this 
critical level is usually too low for the applications of interest. Therefore, it becomes 
essential to explore means for achieving reflectionless propagation of more powerful 
electromagnetic pulses through the atmosphere. Resently, Gildenburg et. al. 8 have point 
out that it is indeed possible to achieve reflectionless propagation of powerful EM pulses in 
an ionizing medium. Their analytical results show that a frequency auto-conversion 
process is in play in a nonstationary plasma. The carrier frequency co of the pulse always 
shifts upward with the growth of plasma frequency cop e under the condition that the 
ionization frequency of the background air is much large than the collision frequency of 
electrons. Therefore, the plasma during breakdown will always remain transparent to the 
pulse (i.e., co > co pe )• Chamber experiments have also been conducted recently to 
demonstrate the frequency auto-conversion phenomenon 9,10 . The experimental results 
show that if the pulse is not distorted seriously (i.e., keeping co > co pe ) at the incident 
boundary, the pulse will not be distorted seriously and will experience frequency up-shift 
while ionizing the air during the subsequent period of propagation through the 
experimental chamber. 

In this paper, a comprehensive theoretical and experimental study of the frequency 
auto-conversion process is reported. Presented in section II is a set of governing equations 
forming the basis of the theoretical model for describing pulse propagation in an ionizing 
air. The equation giving an explicit description of the frequency auto-conversion process is 
also derived and discussed. The pulse propagation experiments and the results are 
presented in section III. A comparison between the experimental results and the theory is 
made in section IV in which a discussion of the work is also given. 

2. PROPAGATION OF A POWERFUL EM PULSE IN AN IONIZING AIR 

The wave equation governing the propagation of an EM pulse in an ionizing 
medium is derived to be : 


2 



( 1 ) 


-c 2 ^-+co pc 2 [l-iv/(co+iv)]}E(z,t) = 0 
3t 2 az 2 


where the electron momentum equation has been used to determine the current density 
induced by the wave field in the plasma; v is the electron-neutral collision frequency, 

v=v+vi , and Vi is the ionization frequency. 

Since the plasma is generated by the pulse, its density (in terms of cope 2 ) is a space- 
time function and is strongly dependent on the intensity and length of the pulse and on 
the loss processes including attachment, recombination and diffusion. When the power 
level of the pulse in consideration far exceeds the breakdown threshold of the air, 
however, the ionization rate vj of the air dominates over all the loss rates of the induced 
plasma. Thus, the rate equation for cope 2 during the pulse period can be simplified to be : 


d 

5-co, 

3 t 


pe 2 = ViCOpe 2 


(2) 


where cope 2 =47cnee 2 / me and the ionization frequency vi is modeled to be 11 

Vi = 3.83 x 10 2 v a [e 3/2 + 3.94 e 1/2 ]exp[-7.546/£] (3) 


in which 8= I E/Eth 1 is the wave field normalized to the breakdown threshold field Eth and 
v a is the attachment rate. 

Equation (1) is analyzed by setting E(z,t)=A(z,t)e 1<!)(z ' t) , where the amplitude function 
A(z,t) varies in space and time much slower than the phase function <{>(z,t). This 
assumption is justifiable for a pulse containing many oscillations. Hence, 
I ainA/at I « I dtydt I and I 31nA/3z I « 1 3<t>/3z I . Using these approximations and the definition 
of local frequency co=9<|)/3t and wavenumber k=-d(p/dz, (1) is reduced to 

3(yA^) + 3(kcW) 2 v ,» (4) 

3t az ^ 

and _ 

CO 2 - k 2 c 2 = COpe 2 [ 1 -V v /(o^+v 2 )] (5) 

Equation (4) is the continuity equation for the energy density of the pulse, while (5) 
is the local dispersion relation. From the definition of co and k, it requires that 

^ + ^ = 0 ( 6 ) 

dz 3t 
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Equations (2)-(6) form a complete set of system equations describing the propagation 
of a one dimensional pulse in a self-generated plasma. To single out the frequency auto- 
conversion process, its governing equation is derived by taking the partial time derivative 
of (5) and using the relation (6) to simplify the result, the resulting equation is 

IF + = ^ 1 " v v/(a 2 +v 2 )]/[ 1 -^vv/f^+v 2 ) 2 ] COpc 2 (7) 


where v g = 9co/3k is the group velocity of the pulse. 

Under the conditions that vv «(co 2 +v 2 ) and cope^f^+v 2 ), (7) is reduced to the one 
derived by Gildenburg et. al. 8 


_d 

dt 



( 8 ) 


where ^ = Jj- + v g^- is the total time derivative along the group trajectory of the pulse. 

We now consider an example which is used to illustrate the physical implication of 
S). Considering a rectangular pulse as shown in Fig. 1, the density of the generated plasma 
{, a given time is expected to grow near exponentially from the background level at the 
; -ading edge of the pulse to a peak value at the tailend of the pulse. Behind the pulse, 
iasma density drops from the peak value at a rate depending on the loss mechanisms. 

he likely density profile is also plotted in Fig. 1. It is shown that dco P e 2 /dz<0 in the pulse 
iterval and, consequently, d(co 2 -co P e 2 )/ dt>0. It, thus, leads to the conclusion that the plasma 
. ill remain transparent (i.e., co pe < co) to the pulse as long as co > co pe at the incident 
oundary. This conclusion seems trivial for a rectangular pulse propagating in a uniform 
medium in which dcope 2 / dt<0. However, it becomes a significant result for pulse 
propagation in a nonuniform medium in which dco P e 2 /dt>0 can occur. An example of this 
case is the propagation of a pulse through the atmosphere. Since the air pressure from the 
ground to the space decreases with the altitude, the breakdown threshold of air also 
decreases with the altitude before reaching an altitude having Paschen breakdown 
minimum. Beyond that altitude, the breakdown threshold increases with the altitude . 
Thus, if the pulse is transmitted from the ground to the space, co P e 2 may increase 
jntinuously along the trajectory of the pulse and exceed <oo 2 =co 2 (t=0,z=0) before reaching 
;e altitude of Paschen minimum. However, there is no cutoff reflection throughout the 
propagation, if co > co P e at the incident boundary. This is because co is also continuously 
-iiifted upward through a frequency upconversion process. One can rewrite (8) to be 
■co 2 /dt = vjco P e 2 >0 , it explains why co also increases along the trajectory of the pulse. It is 
: 1 oted that such a frequency auto-conversion process is also applicable for pulse 
propagation from space to the ground since the altitude of the Paschen breakdown 
minimum is located in between. 
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Fig.l Envelope of a rectangular 
pulse and the pulse-generated 
electron density distribution. 


Microwave 



Fig.2 Experimental setup. 


3. EXPERIMENTS 


The experiment of pulse propagation is conducted in a vacuum chamber filled with 
drv air The chamber is made of a 2 foot cube of Plexiglass. Microwave pulse is fed into the 
cuL through a S-band microwave horn placed at one side of the chamber. A second S- 
band horn placed at the opposite side of the chamber is used to receive the transmitted 
Dulse The spectra of the incident pulse and transmitted pulse are then compared. The 
microwave power is generated by a single magnetron tube (OKH 1448 ) driven by a soft 
tube modulator. The magnetron produces 1MW peak output power at a central frequency 
of 3.27 GHz. The modulator uses a pulse forming network (PFN) having a pulse width of 

us and a repetition rate of 60 HZ. 

The block diagram of the experimental setup is shown in Fig.2. A directional 
coupler connected to the incident horn is used to monitor the spectrum of the input pulse 
which is used as the reference. The air pressure in the chamber is about 0.08 torr, and thus, 
the electron-neutral collision frequency v is below 5xl0*sec-i . Limited by the available 
power of the magnetron tube, the ionization frequency Vi in such a background pressure is 
less than 10»sec-i. Hence, vv«co 2 . We also limit the power level of the input pulse so that 
the plasma density near the incident boundary is below the cutoff density, . • 
cope 2 (z=0)<(co 2 +v 2 ). The assumptions used to reduce (7) to (8) are then justi la e wit 
present experiment. 

The generated plasma density along the path of the pulse is measured in te rms ^o 
the enhanced airglow whose relationship with the actual plasma density is ca 1 X. 

Langmuir double probe. Using a focusing lens to localize ^1.“ ft . 
temporal evolution at a fixed point is then recorded on the oscilloscope thro g 
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photomultiplier tube. By moving the optical probe parallel to the path of the pulse, the 
variation of the peak airglow intensity along the path of the pulse is determined. It is 
noted that the airglow is generated mainly through impact excitation of neutral gas, which 
requires that electron energy exceeds 2eV . Hence, the axial distribution of the peak airglow 
intensity does not really represent, in general, the axial distribution of the peak plasma 
density. Nevertheless, it provides a good approximation especially since only the 
distribution of the peak airglow intensity is involved. Then the absolute value of the 
plasma density can be calibrated by measuring the plasma density with a Langmuir double 
probe at a point where the airglow intensity is also recorded simultaneously. 



CTR 3.2733 GHz SPAN 1 MHz/ RES BW 10'kz VF off 
REF 2.24mv/2.24mv linear for il)/ (2) ATTEN 0 a8 SWP 5 sec/ 


z (Inches) 

(a) Plasma density distribution 
with input power Pi LlP t h. 
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(b) Frequency spectra of the input 
(curve 1) and output (curve 2) pulses. 



(c) The envelope of the input pulse. 


Fig.3 
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We now compare the frequency spectrum of the input and output pulse for two 
different incident power levels, Pi - l.lPih and P 2 - l-46P lh . Although the differences of 
these two levels seem not so significant, the results turn out quite different because the 
ionization rate is extremely sensitive to the field intensity. The corresponding density 
distribution along the propagation path and the envelope of the input pulse for each case 
is also measured for comparison. In the first case. Pi ^ l.lPih the peak density 
distribution of the generated plasma is shown in Fig.3a, which shows that the plasma 
density in this case is peaked at the incident boundary and decays quickly. The second pea 
is probably caused by the reflection of the chamber wall at the exit side. Excluding this pe _ , 
the density distribution appears to have a similar feature as the predicted one presented in 
Fia.l. Correspondingly, the induced frequency shift as shown m Fig. 3b. Presented m F jg. 3c 
is °he envelope of the input pulse. In Fig. 3b, curve 1 is the frequency spectrum of the mput 
pulse while curve 2 is the frequency spectrum of the output pulse. The frequency 
difference between the peaks of these two curves represents the frequency shift attributed 
to the plasma veneration in the chamber. The power level P 2 ~ 1 in the second case 

is chosen for achieving the maximum frequency shift without sacrificing the is or ion 
the pulse The set of frequency spectra for the input pulse (curve I) and output pulse 
Scur^e a'is presented in Fig.4b. As shown in Fig. 4a, the plasma density m this case is 
much higher than that of the previous case shown in Fig. 3a. 



(a) Plasma density distribution (b) Frequency spectra of the input 

with input power P 2 - l-46P lh . (curve 1) and output (curve 2) pulses. 


Fig-4 

4. COMPARISON AND DISCUSSION 
With the aid of (2) and the definition of ^ , (8) can be reduced to 
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(9) 


^<0 2 = ViCOpe 

It is integrated to obtain the frequency shift 

Af = (f 0 /2c)f o L dz Vi(n/nc)/(l-n/n c ) 1/2 (10) 


where dt = dz/v g ~dz/[c(l-n/nc) 1/2 ] and v 2 « coj = co 2 (0) are assumed; L is the length of the 
chamber. 

Therefore, Af can be evaluated once n(z) and vi(z) are specified. Since the length of 
the chamber ~ 0.6m is much shorter than the pulse's spatial length ~300m, the spatial 
variation of vi will be neglected in the following calculation. Moreover, if the input pulse 
is approximated by a rectangular pulse having a constant amplitude Eo and a pulse length 
x, vi=vi(Eo) during the pulse duration is approximately calculated from the relation 
vi=ln(n/n 0 )/t, where n 0 =n(z,t=0 ) is the initial background electron density in the chamber. 

Thus (10) reduces to 


Af ~ [fo /2ct]| q dz (n/n c )ln(n/no)/(l-n/n c ) I/2 (11) 

We now use the density distribution displayed in Figs 3a and 4a to carry out the 
integration of (11). The corresponding results are obtained to be Afi=0.27Mhz and 
Af2=0.40Mhz. From Figs 3b and 4b, the frequency shifts are measured for the two 
corresponding cases to be Afi=0.300Mhz and Af 2 =0.426Mhz, respectively. A good agreement 
between the theoretical and the experimental results is clearly shown. 

In this work, a frequency auto-conversion process for a powerful microwave pulse 
propagating in a self-generated plasma is studied. The theory is first presented and the 
equation explaining the frequency auto-conversion process is derived. This equation 

indicates that if co > co pe at the incident boundary, and the depletion of the pulse energy by 
the self-generated plasma is small, then w increases along the propagation path of the pulse 
and is always larger than co pe throughout the propagation. A pulse propagation experiment 
for confirming the theoretical prediction of frequency auto-conversion phenomenon is 
conducted. The frequency spectra of the input and output pulses are recorded 
simultaneously for comparison. The corresponding peak density distribution along the 
propagation path of the pulse is also determined through the measurement of airglow 
enhancement by an optical probe. Using the measured density distributions, the 
corresponding frequency shifts for two incident pulses of different amplitudes are 
evaluated. They are then compared with the experimental results. A good agreement is 
found. 

The frequency auto-conversion process enables a high power electromagnetic pulse 
to propagate through an inhomogeneous medium, such as the atmosphere, whose 
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breakdown minimum is located between the incidence boundary and the destination. This 
process is thus expected to make the operation of the high power radar and directed energy 
system easier and more effective. 
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Simulation of intense microwave pulse propagation in air breakdown environment 

S.P. Kuo and Y.S. Zhang 
Weber Research Institute 
Polytechnic University 
Farmingdale, NY 11735 

ABSTRACT 

An experiment is conducted to examine the tail erosion phenomenon which occurs to an intense 
microwave pulse propagating in air breakdown environment. In the experiment, a 1 MW microwave 
pulse (1.1a*s) is transmitted through a large plexiglas chamber filled with dry air at about 1~2 ton- 
pressure. Two different degrees of tail erosion caused by two different mechanisms are identified. 
This experimental effort leads to the understanding of the fundamental behavior of tail erosion and 
provides a data base for validating the theoretical model. 

A theoretical model based on two coupled partial differential equations is established to describe 
the propagation of an intense microwave pulse in air breakdown environment. One is derived from 
the Poyntmg theorem, and the other one is the rate equation of electron density. A semi-empirical 
formula of the ionization frequency is adopted for this model. A transformation of these two 
equations to local time frame of reference is introduced so that they can be solved numerically with 
considerably reduced computation time. This model is tested by using it to perform the computer 
simulation of the experiment. The numerical results are shown to agree well with the experimental 
results. 


1. INTRODUCTION 

When a high power microwave pulse propagated through a region of the atmosphere where the 
breakdown threshold field of the background air is lower than the field intensity of the pulse, the 
ionization is expected to occur and produces a space-time dependent plasma. Such a plasma can 
attenuate the pulse and gives rise to a phenomenon called tail erosion which plays the primary role in 
limiting the transmission of the pulse energy. 1 " 4 Moreover, the propagation characteristics of the 
pulse can also be radically modified by the nonlinear and nonlocal effects of the self- generated 
space-time dependent plasma. 4,5 Therefore, any meaningful theoretical effort on the study of pulse 
propagation through the atmosphere requires a self-consistent description of the propagation process 
and a justification with available experimental results. 6 ’. • . 

Basically, there are two fundamental issues to be addressed. The first one concerns the optimum 
pulse characteristics for maximum pulse energy transfer through the atmosphere. The second 
concern is maximizing the ionizations in the plasma trail following the pulse for providing sufficient 
focusing to compensate for wave spreading. In general, these two concerns are interrelated ahd 
must be considered together. This is because the density of the plasma providing the focusing effect 
depends on the intensity of the pulse. On the other hand, the loss of the pulse’s intensity may cause 
excessive ionization in one region where the overdense plasma can cut off the propagation of a large 
portion of the pulse. Once reflection occurs, the remaining pulse will become too narrow for 
effective ionization during the remaining journey. The more serious matter is that the remaining 
pulse will be too narrow to carry enough energy as desired. 

In this research effort, a pulse propagation experiment is first conducted in a chamber to 
demonstrate the tail erosion phenomenon. A theoretical model describing the propagation ot an 
intense microwave pulse in air breakdown environment is then developed. We next use the 
theoretical model to simulate the experiment whose results, in turn, are used to validate the model. 
This model can then be used to simulate the more general cases and determine the conditions 
avoiding the occurrence of cutoff reflection. 


The pulse propagation experiment and results are presented in Section 2. The modal equations 
describing pulse propagation are presented in Section 3. In Section 4, a numerical analysis 
simulating the experiment is made. A frame transformation to convert the modal equations into a 
local time domain description is introduced in Section 5. The resultant equations are shown more 
readily for describing long distance propagation because much less computing time is required. The 
new set of modal equations is normalized according to the relevant parameters of the atmosphere 
and is ready for the numerical analysis of an EM pulse propagating through the atmosphere. This 
work is summarized and discussed in Section 6. 

2. EXPERIMENTS 

Experiments are conducted in a chamber made of plexiglass filled with dry air at about 1 ton- 
pressure. 3 A microwave pulse is fed into the cube by an 5-band microwave horn placed at one side of 
the chamber. A second S-band horn placed at the opposite side of the chamber is used to receive 
the transmitted pulse. The chamber is shielded with microwave absorbers so that the microwave 
reflection from the walls and nearby structures can be minimized. Shown in Fig. 1 is the block 
diagram of the experimental setup. 


screen room 



Figure 1. Experimental setup of high power microwave pulse propagation 


Tail erosion is a common phenomenon appearing in the propagation of a high power microwave 
pulse (HMP). This phenomenon is demonstrated by the snapshots presented in Fig. 2, where 1.1/iS 
pulses, with four consecutively increasing amplitudes, are transmitted through the chamber. The 
first one has a field amplitude below the breakdown threshold field of the background air so nothing 
happens. Consequently, the received pulse shape as shown in Fig. 2(a) is undistorted from that of 
the input pulse. However, once the field amplitude of the input pulse exceeds the breakdown 
threshold, tail erosion occurs and becomes more serious to the larger amplitude pulses, as is shown 
by the subsequent snapshots (Figs. 2(b)-2(d)). This is understandable because the strong increase of 
the ionization rate with field amplitude allows more electrons, which attenuate the pulse, to build up. 
The last two snapshots (Figs. 2(c) and 2(d)), show that the pulses have been strongly eroded. 
However, a clear distinction between the fundamental behavior of tail erosion in the two cases is also 
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noticed In one case corresponding to the third picture (Fig. 2(c)), the erosion to the tail of the pulse 
is not complete. In other words, the received pulse width extends to the original width. In the other 
case (Fis 2(d)) a large portion of the pulse is more or less eroded completely during the finite 
propagation period. Obviously, it is a different mechanism responsible for the second case. The 
ionization frequency becomes so large in the second case that the electron density exceeds the cutoff 
density of the wave before the whole pulse passes through the chamber. The overdense plasma 
screen reflects the rem ainin g portion of the pulse and causes even more severe tail erosion. In order 
to verify that the tail erosion in the second case is indeed caused by cutoff reflection, the reflected 
power tor each case is also measured for comparison. Each of the two snapshots (Figs. 3(b) or 3(d)) 
presented on the RHS of Fig. 3 is the envelope of the reflected pulse correspondmg to each received 
pulse presented on its left (Figs. 3(a) or 3(c)). The second set of snapshots (Figs o(c) and 3(d)) 
shows that strong reflection and complete erosion appear together. In conclusion, two mechanisms 
responsible for the tail erosion are identified. One is due to attenuation by the self-generated 
□erdense plasma. The other one is caused through reflection by the self-generated overdense 

plasma screen. 

Our experimental results (Figs. 3(a) and 3(c)) clearlv demonstrate that an increase of pulse 
amplitude may not help to increase the energy transfer by the pulse. This is because wo tail erosion 
mechanisms have different impacts on causing the degradation of the pulse energy transfer. Ome 
cutoff reflection is allowed to happen, the additional energy added to the pulse is practically wasted 
by the reflection. The practical way to solve this problem is either lowering the amplitude of the 
nulse or narrowing the pulse width so that the excessive ionization can be avoided. It is therefore, 
essential to determine the optimum parameters of the pulse for achieving the most effective energy 
transfer and ionization. Because the effect of pressure gradient and large propagation distance 
cannot be incorporated in the chamber experiment, these parameters can be determined via he 
numerical simulation based on a proposed theoretical model which, however, has to be validated 
first by the available experimental data. 

3. THEORETICAL MODEL 

The Poynting theorem is used to describe the energy flow of the high power wave propagation in 
air breakdown environment, it leads to the equation 


V • — (ExB) + *£- (E 2 + B 2 )/8 tt = -J • E 
4tt ai 

whereT is the wave-induced current density in the self-generated plasma. 

The self-generated plasma is described by the rate equation for the electron density 

= (^i-^a)n - Tfn 2 

at 


(i) 


together with the momentum equation of electrons 


— V = - — - (i/+i/,)V = - — - i/i V 
dt m m 


( 2 ) 


( 3 ) 


where -y is the recombination coefficient, 2 /, v h and y a are collision frequency, ionization frequency, 
and attachment frequency, respectively, and i/ 1 =u+u i . 

In the present study, a one-dimensional model of wave propagation under forward scattering 
approximation is considered. We express the wavefield and the velocity response of the electrons as 

E = line 4 * + c.c. and V = ^ Voe^ 1 + c.c. ( 4 ) 

2 2 

where E 0 and V 0 are the vector phase amplitudes and functions of z and t. 


Substituting (4) into (3), yields 


V 0 = - (z /1 + - iu/)* 1 % 


m 


(5) 


Thus, the induced current density is derived to be 


J = -enV = e 


_ 


ne 2 

2m ( *^ + ¥ )2+ " 2 


En + c.c. 


- e 


-iwt 


^ O'! + -j£ + iw) E 0 + C.C. 


87r(ur 4-r/j) 


( 6 ) 


where we have neglected — terms in the denominator of (6) for a slowly time varying function E 0 , 

d a 

i.e., | ^r/n E o I «w- 


Taking time average of J • E over one period of the wave carrier, yields 

_ y 2 — TT + T aT ' 

87t(o; +z/i) 2 at 


< J-E > = 


(7) 


The time average of (1) then becomes 


d c „ p 5 _1_ 

52 8tt E ° B ° 5t 16 tt 


87t(u; 2 +I/j) 


Bn + 


Eg 


1 + 


w 


,2 > 


u/+v?J 


Eg 


( 8 ) 


where the forward scattering approximation is used to simplify the resultant equation. This 
approximation is used because we are only interested in the energy transfer in the forward direction 
and is justifiable because the density distribution of the induced plasma is, in general, smooth enough 
for neglecting the effect of multiple scattering. 


Introducing the notations 


P = 


1 &K 


Big + 

f 2 > 

l + Wp 

Eg 


w 2 +vfj 



and 


s = i E » B » s i 


_ _1_ 
87T 

, 1/2 


Eg 


(9) 


( 10 ) 


Eg [l-u; 2 /(w 2 +v?)J =v g P 

where P and S are the energy density and power flux of the wave respectively, and 
v g = c[l-cjg/(w 2 +i/ 2 )] 1/2 is the (group) velocity of the energy flow. 

Equation (1) is finally reduced to: 


( 11 ) 



d_ 

dz 


(VgP) = -/TP 


where (3 = ui Wp / (w 2 + u \ ) . 

Equation (11) is coupled with Eq. (2) through the ionization frequency u h which is given by 9 

vj? = 2.5x 10 7 [8.8(P) 1 / 4 + 2.236(P) 3 / 4 ] exp(-7.546/P 1/2 ) (sec-torr)* 1 (12) 

where ? is the background pressure in torr; P=P/P CT , P cr =Ei/87r, and 

E ct =36.84p(l + w 2 /i/ 2 ) 1 / 2 V/cm is the breakdown threshold field of the background air of pressure p. 
These two equations, together with (12) give a self-consistent description of pulse propagation in an 
air breakdown environment. 


4. COMPUTER SIMULATION 


Introducing the dimensionless variahks. and functions: n^n/im P=P/Pcr, t=\/u^u t, 

z' =\/v~v z/c, V g =(l-n) 1/2 , and P=(vi/Vw ) h where n c = (m/4 7 re 2 ) (w 2 +i/i), (2), (11) and (12) 
are normalized to be 


— = (v a M 1/2 An 


( 13 ) 


+ — (v g P) = -/? P 
dr dz % 


(14) 


and 

A = I'i/t'a-l = 3.83 x lO 2 ^ 74 + 3.94P 174 ] expP^b/P 172 ] - 1 (15) 

where the recombination loss term in (2) is neglected and the attachment frequency 
i/ a = 6xl0 4 psec' 1 is measured by the Bragg scattering method. 

We now model the envelope of the normalized energy density of the incident pulses employed in 
the chamber experiments to be 

P 0 = P(r,o) = aexp[-(r/r 0 ) 10 ] (16) 

where t§~\J vv^ to/2 is the normalized pulse width; t 0 is the pulse width and set equal to 1-1 /js in 
this simulation. 


To solve (13) and (14) numerically, we used a software package, developed by Richard F. 
Sincovec and Nile K Madsen 10 for nonlinear partial differential equations. This package is an 
interface routine which uses centered difference approximations to convert one dimensional systems 
of partial differential equations of the form 


d~u 

dt 


= F 


da dV ) 

l ’ x ’ u> aT ’ ax 2 


(17) 


into an explicit system of first order differential equations of the form 

nm = G(t, x,"u ) (18) 

at 

Then it u tiliz es the widely used ODE integrator subroutine 11 to solve these partial differential 
equations by calculating the time derivatives of the functions at mesh points from neighboring spatial 
points and integrating in time domain to get the solution. This program is written in FORTRAN and 
can be run on Cray in the Pittsburgh Super Computing Center. 



Figure 4. Boundary and initial condition of the simulation 

The boundary condition and initial condition for the pulse’s energy density are best shown in Fig. 
4. The initial condition of electron density is n = no for all z; and the density at the boundary will be 
determined from P(t,z = 0) with the restriction of n<l. The differential equations are solved in a two 
dimensional surface determined by z=z' and t=r axes. The number of the mesh points are chosen 
so as to minim ize the numerical error. For the results shown here, we have N m = 1000, N, = 1000. 

Shown in Fig. 5 are four different degrees of tail erosion on four different incident pulses a x = 1, 
a 2 =3.5, 05 =5, °04 = 9 propagating_ over the chamber length (x=66 cm), where the vertical axis 
represents the pulse’s amplitude ( A* =P] / ). Comparing these plots with those presented in Fig. 2 of 
the corresponding experimental results, good agreement is found. 



Figure 5. Tail erosion of 1.1 us microwave pulses for four consecutive 
increasing amplitudes (computer simulation) 
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Figure 6. (a) Microwave pulse envelop as a function of time and space 



Figure 6. (b) Electron density as a function of time and space 
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The num erical results (Fig. 5) of tail erosion phenomenon also reveal two distinctive features 
similar to the experimental observations. For instance, the pulse envelope presented in Fig. 5(b) 
shows a linear decay starting from the front portion of the pulse. In constrast, the pulse envelope 
presented in Fig. 5(d) shows a sharp drop in the amplitude rather than a gradual decay. These 
features confirm with the experimental results shown in Fig. 3, where two received microwave pulses 
with the abovementioned two different features of tail erosion are recorded (Figs. 3(a) and 3(c)). 

The good agreement between the numerical simulation and the experimental results suggests 
that Eqs. (2) and (11) together with (12) represent a very good theoretical model for describing high 
power microwave pulse propagation in air breakdown environment. Thus, the numerical solutions of 
these two modal equations can provide information about the spatial variations of the pulse energy 
and electron density which are not easy to be measured experimentally. 

Shown in Figs. 6(a) and 6(b) are the three dimension drawings of pulse amplitude and electron 
density versus space and time for the case of pulse amplitude a= 15, and the initial electron density 
n o = 10 6 cm' 3 . It shows that the electron density will rise very rapidly at the boundary to reach the 
cutoff density "h= 1, and thus to form a high density thin layer. This high density electron layer will 
absorb and reflect the remaining portion of the pulse. Therefore, only the leading edge of the pulse, 
about one-fifth of the initial pulse width, passes through the high density layer and propagates 
continuously in the forward direction. The transmitted pulse now have lower amplitude and much 
narrower width than those of the incident pulse, the generated electron density, hence, drops sharply 
after the thin layer and maintains a relative constant level following the trail of the pulse. 

5. FRAME TRANSFORMATION OF THE MODAL EQUATIONS 

In validating the theoretical model by the experimental results, one has only to integrate 
numerically the modal equations (2) and (11) over the chamber distance which is relatively short, 
and thus, the computation time does not appear to be a problem. However, in the practical 
applications of pulse propagation over a large distance, the required computation time becomes an 
important issue. Apprarently, the cause of consuming large computation time comes from the 
simultaneous appearance of space and time derivatives in Eq. (11). One common way to eliminate 
such a problem is to use a frame transformation to convert (2) and (11) into a local time domain 
description. In doing so, the space and time derivatives can be separated from the same equation, 
and the computation time on the resultant equations for the same propagation distance can be saved 
tremendously. 

Since the group velocity v g of the pulse is a space and time function, the frame transformation is 
defined to be 

z' ' 

t = t' + / dz'7V g (U,z") 

and (19) 

z = z' 

where t' is the local time of the pulse and V g (t',z') =v g (t,z). Thus, the partial differential operators 
in the two frames are related by 

d/dt = f(t' ,z')d/3t' 

and 

d/dz = d/dz' - [f(t',z')/V g (t',z')]a/at' 
where f(ttyz') = [l-/dz^(V gt '/Vl)r 1 and V*- =dV g /&' -(^Vgity) (dN/dt'). 


Equations (2) and (11) are then transformed to be 

■~"N = [(^.)N^N 2 ]/f 

dt 

and 

^ 7 V g P 1 =-tS-(V gt -/Vg)f]P 1 

az 

where N(t',z')=n(t,z) and Pi(t',z')=P(t,z). 

With the aid of (21) and (22), the integration in f(t',z') can be carried out, and it leads to 

z' 

f(t'.z') = exp{-/dz"c 2 [(i/ 1 -^)N- 1 N 2 ]/2V|n c } 

o 

z" 

= (V g Pi/Vg 0 P 1 „)exp{i/ I / dz'(N/nc)/V g } 

0 

where V g0 =V g (t',z' =0) and Pio =Pi(t',z' =0). Substituting (23) into (21) and (22), yields 

~“N = x- 1 AN(V g0 P I „/VgP 1 )e X p{-t- 1 /dz"(N/n e )/V g } 

ut Q 


( 21 ) 

( 22 ) 


(23) 


(24) 


and 

~Z7 VgP, --[*,+ c 2 J/.A/2^]P,N/n c (25) 

az 

where A is defined the same as (15) witih P =P 1 /P cr ; and z/ x -u+u^. 

Equations (24) and (25) can be solved by using ODE software package directly for pulse 
propagation over a long distance. The required computation time will be much less than that needed 
for solving the original Eqs. (2) and (11) which need an additional interface routine converting (11) 
from the PDE or into an explicit ODE or before using the ODE integrator subroutine. 


The main purpose of this study is to develop a useful model for numerical analysis of pulse 
propagation through the atmosphere. The feature of altitude dependence of the background 
pressure has to be included in the model. This, in turn, causes u, i/ a , u u and all become 
altitude dependent. Therefore, (24) and (25) cannot be normalized directly in terms of these 
parameters as made in Section 4. Instead, the parameter values at the location of 1 torr pressure 
(near the Paschen minimum) are chosen for the normalization purpose. We assume that the air 
pressure at 50 km height is 1 torr, and the pressure is an exponential decay function of the altitude 
z'. Thus, p(z')=760e‘ 1 ' 33>< ^~* z torr, where z' is measured in meters. Let u 0 =u(z' =0), then 
ixlO^z j n terms 0 f the ground level dc breakdown threshold power P co =P <T iV =0,w=0), 
f - • --2.66x10 Z ^ u/hpi-A hP 1 ^irm'liirlv n _ = n _ _ (n ■+• #»'*' ^0 z ^ where 




- 1.33 > 


then P cr (z') = P c0 (77+e' 


), where r)=uf Similarly, n c =n co (?j+e' 


n c0 =nm 2 /47re 2 . At 1 torr pressure, =i/ ac = 6xl0 4 sec' 1 and i/=i/ c = 6xl0 9 sec' 1 for 1 eV electron 
temperature. Including the effect of electron energy change by the incident pulse and assuming that 
the unperturbed background temperature is about 0.01 eV, a modal function for v is given by 
z/=z/ c p[P 1/2 + 0.01] 1/2 , where T=P/P cr = (P/P co ) (rj+e’ 2 - 66 * 10r>z ' Y 1 . We now introduce the 
dimensionless variables and functions as follows: WV->t, (i/ c p ac ) 1/2 z'/c- + z, 

— ' T/ /XT '- N1/ ' - -- J — - h(z) = (^+e' 4 ' 2x W 2 )‘ 1 > Pi=P/P C0 -Pi. 


n = N/nc = (N/nco )h(z)^n. 


where 




Q = (l-n) 1/ ' 2 P 1 -»-(l-n) 1 / 2 P 1 , Pi =Q/(l-n) 1 / 2 , and P-+Pih(z) = (Q/V 1-n )h(z). In terms of these 
dimensionless variables and functions, and notice that (VcA'ac) 1 ^ =316 and (^ac) 1/2 = 1.9xl0 7 sec' 1 . 
Equations (24), (25) and (15) become 

p. z 

— n = 2.41g(z)An(Q o /Q)exp{-2.4xl0 5 ff(x)dx} (26) 

ot 


-^Q = -2.41g(z)j (316) 2 [(P 1 h)V 2 +0.01] + 

and 

A = S^lO^^h) 3 / 4 + 3.94(P 1 h) 1/4 >exp{-7.546/(P 1 h) I / 2 >-l (28) 

where gCz)*^ 2 ' 1 " 10 ” 1 , f(x)=g(x){[P 1 (x)h(x)]‘/ 2 + 0.01]V 2 [n(x)/Vl-n(x)] and Q„=Q(t,z=0). 

Equations (26)-(28) are the set of modal equations describing microwave pulse propagation 
through the atmosphere either from ground to space or from space to ground. If the pulse is oblique 
propagation with an an gle 6 with resp ect to the vertical axis, what one has to do is to replace V g by 
V o cos0, i.e., replacing vk by v 1-n cos# in (26) and (27), for a correct description. In terms of z 
and t, the red distance x and time r are given by x = 15.8 z(m) and r =52.63 t(ns). 

6. SUMMARY AND DISCUSSION 

The propagation of a high power microwave pulse through a self-generated plasma is studied. A 
chamber experiment is first conducted to explore the unique features appearing in the pulse as the 
result of the breakdown of the background air. The most apparent one is tail erosion. It 
demonstrates that air breakdown can indeed radically modify wave propagation and limit the energy 
transfer by the pulse. Two mechanisms responsible for two different degrees of tail erosion have 
been identified. The serious one happens when the initial intensity of the pulse is too high and 
causes complete erosion to a large tail portion of the pulse during the finite propagation. It is 
identified as a result of the reflection by the self-generated ovqrdense plasma screen. It can be 
realized by the fact that the ionization frequency becomes too large, the electron density can build up 
quickly and exceed the cutoff density of the wave before the whole pulse passes through. . The results 
show that once reflection occurs, the remaining pulse will become too narrow to deliver enough 
energy as desired. Therefore, there is a need to find the optimum pulse parameters for achieving 
maximum energy transfer. Consequently, a theoretical model developed for this purpose would be 
more practical. The experimental results of the present work are then used for validating the 
established theoretical model of high power microwave pulse propagation in the atmosphere. 

The self-consistent description of the propagation process is provided by the set of two modal 
equations including Poynting’s equation for the continuity of the power flux of the pulse and the rate 
equation of the electron density. A foiward wave approximation model is used to simplify Poynting’s 
equation and a semi-emperical formula for as a function of the field amplitude is used to express 
the electron rate equation explicitly. This set of modal equations is then used to simulate the 
chamber experiment. The numerical results are shown to agree very wed with the experimental 
results, especially the unique features of two distinctive degrees of tad erosion identified from the 
experimental results can be reproduced by the numerical results of the modal equations. This 
comparison serves as a validation of the established theoretical model. 



In order to improve the numerical efficiency of the model in terms of the required computation 
time and available subroutines for numerical analysis of pulse propagation over a long distance, a 
transformation to the frame of local time of the pulse is introduced. The effect of mace-tune 
dependence of the group velocity of the pulse is included in this properly designed transformation. 
The resultant equations are reduced to the forms which can be solved directly by the available 
subroutine of ODE solver. In terms of relevant parameters of the atmosphere, these two equations 
are normalized for numerical analysis of pulse propagation m the atmosphere It is noted that the 
inhomogeneous feature of the background pressure is preserved m the model, thus, the model is 
capable for a description of long distance pulse propagation. 
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